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ABSTRACT 

We present a review of our recent results for the comparative evaluation of the induced exchange interaction 
and quantum noise mediated by the bosonic environment in two-qubit systems. We report new calculations 
for P-donor-clcctron spins in Si-Ge type materials. Challenges and open problems arc discussed. 
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1. INTRODUCTION 

In this review we present results of our recent investigations 1,2 addressing several problems in the field of 
open quantum systems; some have a rather long history. With the actual experimental probes now being 
carried out at the nanoscale, these problems have become more pressing, and some have actually been 
suggested by the experimental developments. 3-14 Perhaps the most fundamental (and difficult) of these 
problems is the matter of accounting within a calculationally tractable approach for relaxation vs. coherent 
dynamics in open quantum systems. In particular the coherent and quantum-noise features induced by 
the bosonic environment (bath of noninteracting bosonic modes) are especially interesting, because this 
model is widely applicable for quantum computing systems. Here we consider bath-induced RKKY-type 
exchange interactions in two-qubit systems. 

This paper is organized as follows. In the next section, Sec. 2, we formulate the model and discuss the 
issue of the initial conditions. An example of the semiconductor (Si-Ge type) based qubit model is given 
and the nature of the interactions is discussed. 

Section 3 is devoted to the initial dynamics of the induced interaction and quantum noise. In particular, 
in Sec. 3.1 we present an exactly solvable model to obtain the exchange interaction as well as the time- 
dependent correction due to the initial conditions. The model also demonstrates the delay in the response 
of the qubits caused by the finite propagation velocities of the mediating virtual bosons. In Sec. 3.2, 
we study the interplay of the coherent dynamics and quantum noise by calculating the concurrence. 15, 16 
Evolution of the reduced density matrix elements is also investigated. 

At large times, as the system forgets the initial state, the exchange interaction becomes stationary. This 
regime is discussed in Sec. 4. We utilize the master-equation formulation of Markovian dynamics, Sec. 4.1, 
and go over an illustrative one-dimensional (ID) example of the emerging interaction and quantum noise in 
Sec. 4.2. Generalization to higher dimensionality is given in Sec. 4.3. A solid-state based three-dimensional 
(3D) example of a semiconductor qubit model is presented in Sec. 4.4. In the latter discussion we compare 
bath-induced coupling and noise to other interactions, such as an electromagnetic coupling of P-impurity 
spins (qubits). Finally, brief summary and outline of some open challenges arc offered in Sec. 5. 
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2. THE MODEL AND INITIAL CONDITIONS 



We consider two qubits, i.e., two two-state quantum systems, modeled by localized electron spins (1/2) 
immersed in a common quantum environment. The qubits are located at the distance d from each other, 
far enough so that the direct overlap of the electron wave functions is negligible. The environment is 
modeled by a bath of bosonic modes, which are maintained at temperature T. It is usually assumed that 
external influences, as well as possibly internal bath-mode interactions, set a fast time scale t B for the bath 
correlations to equilibrate. The bath modes are then regarded as otherwise noninteracting. At least for the 
low-frequency bath modes, it is usually argued that such a thermalization time for a generic case should 
be of order h/k^T. The thermal-state density matrix of the bath modes, taken here as noninteracting 
bosonic fields, is 



PB 



cxp 



-y^k,g«k, f Qk,;/fcBr 



(1) 



e ^k/k B Ty j s linearly coupled 

(2) 

where the superscripts (j = 1,2) label the two spins, and the bath operators arc given by 



where we set H = 1, and the partition function is Z = 1/ Ilfc(l 
to each qubit, 

r sb = E E <xL > 

j— 1,2 m—x,y,z 



x 3 



The overall Hamiltonian of the qubit-bath system is 



H = H S + H B + H 



SB 



where 



H B = E a;k '£ a k,g ak .£ 

k,C 



and the two spin-1/2 (qubits) will be assumed split by external magnetic field, 

H s = A(al+<j 2 z )/2. 



(3) 



(4) 



(5) 



(6) 



Here A is the energy gap between the up and down states for spins 1 and 2. A natural example of such 
a system are spins of two localized electrons interacting via lattice vibrations (phonons) by means of the 
spin-orbit interaction. 17-20 Another example is provided by atoms or ions in a cavity, used as two-state 
systems interacting with photons. 21 

Our emphasis here will be on calculating and comparing the relative importance of the coherent (induced 
interaction) vs. quantum-noise effects of a given bosonic bath in the two-qubit dynamics. We do not include 
other possible two-qubit interactions in such comparative calculation of dynamical quantities. In Sec. 4.4 
we also include the direct electromagnetic (EM) coupling for a comparison of the bath-induced and EM 
dipole-dipole interaction strengths. 

An interesting realization of the model formulated above can be found in the impurity electron spin 
dynamics in semiconductors. Let us consider P donor impurities embedded at controlled positions in an 
otherwise very clean Si (or Ge) crystal matrix. The system is maintained at very low temperature, as 
appropriate for quantum computing. Therefore, the outer donor-impurity electron remains bound. The 
spins of such localized electrons can be utilized as qubits. They are subject to external magnetic field H, 



which produces the Zeeman splitting (6). The spin-orbit interaction couples the spin to the deformation 
potential fluctuations of the host semiconductor, producing the energy change 



Hi , (7) 

m,Z — x,y,z 

where is the Bohr magneton, and H = {H x , H y , H z }. Here the tensor g m i is sensitive to lattice 
deformations. It was shown 19 that for the donor state which has tetrahcdral symmetry (which is the case 
for P in Si or Ge), the Hamiltonian (7) yields the spin-deformation interaction of the form 

Hi = A^b [e xx <r x H x + e yy a y H y + e zz a z H z + A(a ■ H)/3] + B^ B [e xy {o x H y + a y H x ) + c.p.] . (8) 

Here c.p. denotes cyclic permutations and A is the effective dilatation. The tensor already includes 
averaging of the strain with the gradient of the potential over the donor ground state wave function. 

As before, let us assume that we have two impurities separated by distance d. For definiteness, one can 
direct the magnetic field along the z-axis. Then the spin-deformation interaction Hamiltonian simplifies 
to 

Hi = k^ B e zz a\H z + B/ib {e yz a l y H z + e zx a l x H z ) . (9) 
In terms of the quantized phonon field, we have 17 ' 22 



^ = 2 

where in the spherical donor ground state approximation 18,22 



- /(k) V 8pV^T ^ k ' ikj + ^ ,jk ^ (°* ,£ + ak ' ? ) ' (10) 



/(k) = n/n^ 2 ' 

(l + a 2 fc 2 ) 

Here Ob is half the effective Bohr radius of the donor ground state wave function. In an actual Si or Gc 
crystal, donor states are more complicated and include corrections due to the symmetry of the crystal 
matrix including the fast Bloch-function oscillations. However, the wave function of the donor electrons in 
our case is spread over several atomic dimensions (sec below). Therefore, it suffices to consider "envelope" 
quantities. Thus, in the spin-phonon Hamiltonian (2)- (3) coupling constants will be taken in the form 



9™,t = ~, ~T: on2 (£k,z fc m + £k,mfe) , ( 12 ) 

(1 + a 2 k 2 ) V 8 P Vuj k,a 

where D x = D y = Bfi B Hz and D z = A/ibH z . 

It will be instructive to consider a onc-dimensional calculation, which simplifies the notation. General 
results 1 ' 2 are presented later in this article. One can think of a ID channel geometry along the z direction. 
This will give an example of an Ohmic bath model discussed later. In a ID channel the boundaries 23 can 
approximately quantize the spectrum of the phonons along x and ?/, depicting the density of states except at 
certain resonant values. Therefore, the low-frequency effects, including the induced coupling and quantum 
noise, will become effectively one-dimensional, especially if the effective gap due to the confinement is 
of the order of to c . The frequency cutoff comes from (11), namely, it is due to the bound electron wave 
function localization. A channel of width comparable to ~ a B will be required. This, however, may be 
difficult to achieve in Si or Ge with the present-day technology. Other systems may offer more immediately 
available ID geometries for testing similar theories, for instance, carbon nanotubes, chains of ionized atoms 
suspended in ion traps, 24-26 etc. In the ID case, the longitudinal acoustic (LA, ||) phonons will account 



for the <7£ || component of the coupling, whereas the transverse acoustic (TA, _L) phonons will affect only 
the x and y spin projections. 

It will be shown later that the contributions of the cross-products of coupling constants, g^^g^e)* 
with m 7^ m', to quantities of interest vanish. The combinations that enter the dynamics are 

I z i2 _ AV|ff 2 fc 2 2 _ y 2 _ B 2 p 2 i H z fc 2 

4 P Kw MI (l + a|fc|) 4pKw fci _L (l + a |fc2) 

With the usual assumption for the low-frequency dispersion relations w& i| ps C||fc z and w/c^x ps c±k z , the 
expressions (13) lead to the Ohmic bath model. The shape of the frequency cutoff resulting from (13) is 
not exponential. However, to estimate the magnitude of the interaction one can equivalcntly consider the 
exponential cutoff, 27 with the summation over bosonic modes carried out as 

oo oc 

E "> / <M0 m M| 2 TM = Jdcj a>"cxp(- W K) , (14) 

k <? 

where T(cu) is the density of states. The coupling constants should then be taken as 

z _ a Vb#z „,x y _ B f^ H z 



(15) 



where 5 is the cross section of the channel, and the cutoff is u> c — » ch/ob for the z component, and 
u c — > cj^/aB for the x and y components. 

For numerical estimates, we note that a typical value 18, 19, 22 of the effective Bohr radius in Si for 
the P-donor-electron ground state wave function is 2ae = 2.0 nm. The crystal lattice density is p = 
2.3 x 10 3 kg/m 3 , and the g-factor is g* = 1.98. For an order-of-magnitudc estimate, we take a typical 
value of the phonon group velocity, c s = 0.93 x 10 m/s. The spin-orbit coupling constants in Si are 18, 19, 22 
A 2 « 10 2 and B 2 ps 10 _1 . In the Ge lattice, the spin-orbit coupling is dominated by the non-diagonal 
terms, 18,19 ' 22 A 2 ps and B 2 ps 10 6 . The other parameters are 2a B = 4.0 nm, p = 5.3 x 10 3 kg/m 3 , 
c s = 5.37 x 10 3 m/s, and g* = 1.56. This results in a much stronger transverse component interaction. In 
both cases the magnetic field will be taken of order H z = 3 x 10 4 G. One could use other experimentally 
suggested values for the parameters, such as, for instance, ob- This will not affect the results significantly. 

As mentioned earlier, a realistic phonon environment includes the time scale at which it is reset to the 
thermal state. One can think of the short- and long-time dynamics in reference to this time scale. In the 
long-time regime, it is customary to introduce Markovian-type assumptions, which include resetting the 
bath to the thermal state instantaneously. We will use this approach in the master-equation formulation 
of the large time dynamics, Sec. 4.1. 

For the short-time dynamics, one has to address the matter of the initial conditions at t = 0. We 
will use the initially factorized density matrix, with the thermal state assumed for the bath. In quantum 
computing applications, such an initial condition has been widely used for the qubit-bath system, 28-30 

p(0) = ps(0)p B ■ (16) 

This choice allows comparison with the Markovian results, and is usually needed in order to make the 
short-time approximation schemes tractable; 30 specifically, it is necessary for the exact solvability of the 
model considered in the next section. 

A somewhat more "physical" excuse for the factorized initial conditions is formulated as follows. Quan- 
tum computation is carried out over a sequence of time intervals during which various operations are 



performed on individual qubits and on pairs of qubits. These operations include control gates, and mea- 
surements for error correction. It is usually assumed that these "control" functions, involving rather strong 
interactions with external objects, as compared to interactions with sources of quantum noise, erase the 
fragile entanglement with the bath modes that qubits can develop before those time intervals when they 
are "left alone" to evolve under their internal (and bath induced) interactions. Thus, for evaluating relative 
importance of the quantum noise effects on the internal (and bath induced) qubit dynamics, which is our 
goal here, we can assume that the state of the qubit-bath system is "reset" to uncorrelated at t = 0. 



Let us consider relatively short time scales and analyze how the exchange interaction, accompanied by 
the quantum noise, sets in. One can argue that such an investigation is only feasible provided one knows 
the initial condition for the qubits-bath density matrix. Indeed, the short-time interaction (and quantum 
noise) should depend on the initial condition. The factorized initial condition, with the thermal-state bath, 
assumed here, is expected in most quantum computing applications, as discussed above. 

3.1. Development of Induced Exchange Interaction 

When the evolution of isolated qubits is slow with respect to the other time scales such as that of deco- 
herence, so that one can assume vanishing qubit splitting energy, A = 0, and if only one system operator 
enters the qubit-bath interaction, then one obtains an exactly solvable model. A more general "adiabatic" 
system that allows exact solvability is obtained when Hs commutes with Hsb- For example, for electron 
spins of P-impuritics in Si one has the spin-phonon interaction dominated by such an adiabatic term, and 
only one such system operator in the interaction matters for the onset dynamics. Though this situation is 
rarely the case in quantum computing systems, the present model provides a convenient tool for evaluat- 
ing the initial dynamics of the exchange interaction build-up, as well as for analyzing the response delay 
due to the finite speed of the mediating virtual bosonic particles. These features arc difficult to capture 
analytically in other models. 

Therefore, let us presently take ot v n = ajj = 0, while ^ 0. With the above assumptions, one can 
utilize the bosonic operator techniques 31 to obtain the reduced density matrix for the system (2)-(3), 



where the projection operator is defined as P\ = | Ai A 2 ) (A1A2I, and \Xj) are the eigenvectors of a J x . The 
exponent in (17) consists of the real part, which leads to decay of off-diagonal density-matrix elements 
resulting in decoherence, 



3. ONSET OF CORRELATIONS BETWEEN QUBITS 




(17) 



Re£ AA ,(i) = -^G fe (i,T) (A' 1 -A 1 ) 2 + (A^-A 2 ) 2 + 2 cos 



k 



( 



w fc |d| 



) 



(Ai - AO (A' 2 - A 2 ) 



(18) 



and the imaginary part, which describes the coherent evolution, 




(19) 



Here we defined the standard spectral functions 



27,28 




(20) 



Figure 1. The magnitude of the time-dependent Hamiltonian corresponding to the initial correction as a function 
of time and distance. The Ohmic (n = 1) case is shown. The inset demonstrates the onset of the cross-qubit 
interaction on the same time scale. 



and 



C k (t) = 2^- (oj k t - sinwfci) 



(21) 



Calculating the sums by converting them to integrals over the bath-mode frequencies us in (18) and (19), 
assuming the Ohmic bath with n = 1, for T > one obtains a linear in time, t, large-time behavior for 
both the temperature-dependent real part and for the imaginary part. The coefficient for the former is 
~kT, whereas for the latter it is ^us c . For the super-Ohmic models, n > 1, the real part grows slower, as 
was also noted in the literature. 28, 32-34 

First, let us analyze the coherent evolution part in (17), namely, the effect that the imaginary part of 
£aa' (t) has on the evolution of the reduced density matrix, since this contribution leads to the induced 
interaction. If we omitted (18) from (17), (19), and (21), we would obtain the (coherent) evolution operator 
in the form 

e -i{H int +F{t))t _ (22) 

The interaction H lnt comes from the first term in (21), 



2<r(n)c> c " 

z/2 



n arctan 



(23) 



This expression gives the constant interaction that is important in large time dynamics. We will obtain 
this induced exchange interaction latter using different techniques for more general cases, see Sec. 4. The 
operator F(t) is given by 



, D(us)\q(ui)\ sinuit 
dus y — n cos 



ust 



(24) 



It commutes with H lnt (and with itself at different times), and therefore d (tF(i)) /dt can be viewed as the 
initial time-dependent correction to the interaction. In fact, this term controls the onset of the induced 
coherent interaction; note that F(0) = —H int , but for large times F(t) ~ a*u>™ / (us c t) n . 



Let us consider in detail the time dependent correction Hp{t) = d (tF(t)) jdt to the interaction Hamil- 
tonian during the initial evolution, 



H F (t) 



alala n T(n) u(u c \d\/c s 



bj c t) + u(uj c \d\/c s + tu c t) , 



(25) 



where u(£) = cos[n arctan(^)]/(H-^ 2 ) n//2 . The above expression is a superposition of two waves propagating 
in opposite directions. In the Ohmic case, n = 1, the shape of the wave is simply = 1/(1 + £ 2 ). 

In Figure 1, we present the amplitude of Hp(t), defined via Hp(t) = Hpala^,, as wen as the sum of 
Hi nt and Hp (t) , for n = 1. One can observe that the "onset wave" of considerable amplitude and of shape 
u(£) propagates once between the qubits, "switching on" the interaction. It does not affect the qubits once 
the interaction has set in. One can also see in Figure 1, as well as from (25), that the interaction between 
the qubits is delayed by the time ~ |d|/c s , which is required for the mediating (virtual) bosons to carry 
the response to the other qubit. At the same time, relatively slow decay of the initial correction F(t) may 
necessitate a discussion regarding the meaning of the "coherent" large-time induced interaction (when the 
noise effects have also fully set in). We will offer additional comments in the concluding section. 

3.2. Initial Dynamics of the Density Matrix Elements and Entanglement 

Let us now take the entire solution (17) which includes both the induced interaction discussed above and 
the noise (18). In the exact solution of the short-time model, the bath is assumed to be thermalized only 
initially. At large enough times the externally induced equilibration of the bath should be considered. We 
will account for this later when discussing the pcrturbativc Markovian approach. Here we continue our 
investigation of the short-time model which assumes that the two qubits coupled with the (nonintcracting 
with each other) bath modes constitute a closed quantum system. 

To quantify the interplay of the effects of the induced interaction vs. noise, we evaluate the con- 
currence, 15, 16 which measures the entanglement of the spin system and is monotonically related to the 
entanglement of formation. 35, 36 For a mixed state of two qubits, ps, we first define the spin-flipped state, 
Ps = <J y <J y P*s ^y^y' ano - then the Hermitian operator R = yJ~^J = pWps\/p = s \ with eigenvalues The 
concurrence is then given 16 by 



Considering the Si channel geometry for the interaction, introduced above as an illustrative example, we 
arrive at an approximately adiabatic Hamiltonian (a \ 3> c^' a ) with Ohmic- type coupling. The dynamics 
of the concurrence is presented in Figure 2, and we note that the peak entanglement can reach a sizable 
fraction of 1. The coupling constant a\ is quite small due to the weakness of the spin-orbit coupling of 
P-impurity electrons in Si, which results in low magnitude of the induced interaction (and the noise due 
to the same environment). Nevertheless, one still observes decaying periodic oscillations of entanglement, 
which indicate that an approximately coherent dynamics can develop, due to the bath-induced interaction, 
for several dynamical cycles. 

To understand the dynamics of the qubit system and its entanglement, let us continue with the analysis 
of the coherent part in (17). For each spin, we define the two states |±) = [ ||) ± | J,)] /y/2. After the interac- 
tion, H ln t, sets in (note the time scales in Figures 1 and 2), it will split the system energies into two degen- 
erate pairs Eq — E\ = —Hint and E2 = E3 = H m t- The wave function is then \ip(t)} = exp(— iH ni tt)\i^(0)) . 

For the initial "++" state, |^>(0)) = |++), it develops as \tp(t)} = |++) cosHi„tf + | ) isinHmtt, where 

Hmt = ^intf'm' 7 ™- One can easily check that at times tp = 7r/47ii nt , 37r/47ii nt , ■ ■ maximally entangled 
(Bell) states are obtained, while at times to = 0, 7r/27ij nt , ir/Hint, ■ ■ ■, the entanglement vanishes; these 
special times can also be seen in Figure 2. 
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(26) 



Figure 2. Development of the concurrence as a function of time, calculated with af = 0.5 • 1CP 7 and hsT/cuc — 
0.5 ■ 10~ 2 , which corresponds to the magnetic field H, — 0.53 T and temperatures T = 0.34 K. The left inset shows 
the distribution of the concurrence in the |d|-t plane. The right inset presents the dynamics of the diagonal density 
matrix elements P++ = {+ + |/£>s(i)| + +), etc., on the same time scale. 



The coherent dynamics just described is only approximate, because the bath also induces decoherence 
that enters via (18). The result for the entanglement is that the decaying envelope function is superimposed 
on the coherent oscillations described above. The magnitudes of the first and subsequent peaks of the 
concurrence are determined only by this function. As temperature increases, the envelope decays faster 
resulting in lower values of the concurrence. 

Note also that the non-monotonic behavior of the entanglement is possible only provided that the initial 
state is not trivial with respect to the induced interaction, otherwise only the phase factor is developed. 
For example, taking the initial state [ | — h) + |H — )] /V%, in our case would only lead to the destruction of 
entanglement, i.e., to a monotonically decreasing concurrence, similar to results of other studies. ' 

For the model that allows the exact solution, i.e., for Hs = 0, one can notice that there is no relaxation 
by energy transfer between the system and bath. The exponentials in (17), with (18), suppress only the off- 
diagonal matrix elements, i.e., those with A 7^ A'. It happens, however, that at large times the d-dependence 
(the cosine term in Re£AA') is not important in (18). Therefore, one can show that Re£AA'(i — * 00) 
vanishes for certain combinations of A 7^ A'. Specifically, the limiting t — > 00 density matrix for our initial 
state (|H — h)) retains some non-diagonal elements, 
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1 1 
1 1 

3 / 



p{i — > 00) 








(27) 



V -1 



The basis states here are |++), \-\ — ), | — h), and | ). The significance of this and similar 39 results is 

in the fact that in the model with Hs = and non-rethermalizing bath not all the off-diagonal matrix 
elements need be suppressed by decoherence, even though the concurrence of (27) is zero. 



4. INTERPLAY BETWEEN THE INDUCED INTERACTION AND QUANTUM 

NOISE AT LARGER TIMES 



4.1. Master-Equation Approach 

In this section we present the expressions for the induced interaction and also for the noise effects due to 
the bosonic environment, calculated pcrturbatively to the second order in the spin-boson interaction, and 
with the assumption that the environment is constantly reset to thermal. 2 

The dynamics of the system can be described by the equation for the density matrix, 

ip(t) = [H,p(t)]. (28) 

In order to trace over the bath variables, we carry out the second-order perturbative expansion. This 
dynamical description is supplemented by the set of Markovian assumptions, 27, 31 ' 33 ' 40 one of which invokes 
resetting the bath to thermal equilibrium, at temperature T, after each infinitesimal time step, as well as 
at time t = 0, see (16), thereby decoupling the qubit system from the environment. 27,33 This is a physical 
assumption appropriate for all but the shortest time scales of the system dynamics. 28-30 It can also be 
viewed as a means to phenomenologically account in part for the randomization of the bath modes due to 
their interactions with each other (anharmonicity) in real systems. This leads to the master equation for 
the reduced density matrix of the qubits, ps(t) = Tr B p(t), 

oo 

ip s (t) = [H s ,p s (t)] - i J dt'Tr B [H SB , [H SB (t' - t),p B p s (t)]] , (29) 
o 

where H sb (t) = e i{ - HB+Hs "> T H SB e^ HB+Hs ^ T . Analyzing the structure of the integrand in (29), after 
lengthy calculations 2 one can obtain the equation with explicitly separated coherent and noise contribu- 
tions, 

ips (t) = [fleff, Ps (*)] + iMps (t) • (30) 
The effective coherent Hamiltonian iJ ff is 

ffcfF H S +J2 2xT(d)alal - x * (d) (^a 2 + o^) + X «(d) (^a 2 + a 2 ^) - [rf a (0) + ^Jf(O)] {a\ + a 2 ) . 

m— x,y,z 

(31) 

The expressions for the amplitudes \T (d) > XT (d) > V™ (d) , and rj™ (d) will be given bellow. The first three 
terms following H$ constitute the interaction between the two spins. We will argue below that the leading 
induced exchange interaction is given by the first added term, proportional to xT (d) ■ The last term gives 
the qubit Lamb shifts. 

The otherwise cumbersome expression for the noise terms can be represented concisely by introducing 
the noise superopcrator M, which involves singlc-qubit contributions, which are usually dominant, as well 
as two-qubits terms, 

M = £[M4(0)+^M«(d)], (32) 

where the summations are over the components, m = x,y, z, and the qubits, i, j = 1,2. One can define the 
superoperators L a (Oi)0 2 = {O x , 2 }, L(O u 2 )0 3 = 1 O s 2 , and L ± {O u 2 ) = L{0 1 ,0 2 )±L{0 2 ,O x ), 
to write 

M^(d) = C(d) [2Z(c4,<)-Z a (cr>4)] +r? s ™(d) [l+(<4, - L a {a l m ^)\ - i x ™(d)L_((4, , 

(33) 



where we denote qj, = | [cr| , cr^] , and 



(34) 



%{o) = C(o) [2i(< )C 74) +c(o)i+«,^)-^r(o) [l_(<,o + l (< ? 4) 

The amplitudes in (31), (33), and (34), calculated for the interaction defined in (2)-(3), are 

oo 

v"Vrn- V [ \n m I 2 ^COs(k-d) 

Xc ( } ~ ^JT^f W c 2 -A 2 (1-<W) ' (35) 



< l ( d ) = ? E / 7^5 ~th ^ cos (k • d) E * K>« + C 1 " <W)*A) > (36) 

2 Y •/ (2tt) 2fc B T ^1 

> — OO 4 

OO 

xr(d) = -(l-<W)E / ^|< c | 2 cos(k.d) £ ^K,« + gA), (37) 

OO 

,„,,> , ,0 / Vd]i I ml 2 „ ^k,£ Acos(k-d) 

^ (d) = (1 - W J y ^ fe| coth^-i, A2 . (38) 

Here the principal values of integrals are assumed. 

Note that x™(d) appears only in the induced interaction Hamiltonian in (31), whereas f?™(d), x™(d), 
and f?™(d) enter both the interaction and noise terms. Therefore, in order to establish that the induced 
interaction can be significant for some time scales, we have to demonstrate that x™(d) can have a much 
larger magnitude than the maximum of the magnitudes of 77™ (d), x™(d), and ?y™(d). The third and 
fourth terms in the expression for the interaction (31) are comparable to the noise and therefore have no 
significant contribution to the coherent dynamics. 

4.2. Qubit Coupling in a ID Channel 

It is instructive to start with the simple ID example, leaving the derivations for higher dimensions to 
the next subsection. The ID geometry can actually be natural for certain ion-trap quantum-computing 
schemes, in which ions in a chain are subject to Coulomb interaction, developing a variety of phonon-mode 
lattice vibrations. 24-26 

In our case we allow the phonons to propagate in a single direction, along d, so that k ■ d = k |d|. 
Here, for definitcness, we also assume the linear dispersion, ujk = c s k. Furthermore, we will ignore the 
polarization index in the coupling constants: 5^ — > g m (tu). 

The induced interaction and noise terms depend on the amplitudes x™(d), 77™ (d), x™(d), and 77™ (d), 
two of which can be evaluated explicitly for the ID case, because of the 5- functions in (36), (37). However, 
to derive an explicit expression for X™(d) and 77™ (d), one needs to specify the w-dependence in |g m (a;)| . 
For the sake of simplicity, in this section we approximate |c/ m (u;)| 2 by a linear function with superimposed 
exponential cutoff. For a constant ID density of states V/2ttc s , this is a variant of the Ohmic-dissipation 27 
model, i.e., (14) with n = 1 (the case when n > 1 is considered in the next subsection). 

In most practical applications, we expect that A <C w c . With this assumption, we obtain 

X?(d) = ^ (39) 

1 + («c |d|/c s ) 2 



x ™(d)=a5> c £(l-* TOl ,) — cos^M, (40) 
2 w c c s 

C(d) = <^ (1 - S m , z ) - coth -A- cos A_I . (41) 
2 w c 2fc B i c s 

The expression for rj s (d) could not be obtained in closed form. However, numerical estimates suggest 

that 77s (d) is comparable to 77™ (d). At short spin separations d, 773(d) is approximately bounded by 

— a™ A In jj- exp(— ^r^), while at lager distances it may be approximated by a" sin coth 2fc ^ r ■ 

The level of noise may be estimated by considering the quantity 

M = max|C(d),xr(d),C(d)| • (42) 

|d| 

The interaction Hamiltonian takes the form 

ffint = - - - 2 — ~ 2 £ a l^W m . (43) 

1 + (u> c \d\/c a ) m ^ y>z 

This induced interaction is temperature independent. It is long-range and decays as a power law for large 
d. Note that the expression (43) is the same as the induced interaction obtained from the short-time model 
(23), provided the proper spin-phonon interaction component is considered. 

If the noise term were not present, the spin system would be governed by the Hamiltonian H$ + -Hint- 
To be specific, let us analyze the spectrum, for instance, for af = ot\, a* = 0. The two-qubit states would 
consist of the singlet OtlHIt/O/V 2 " and the split triplet |Tt), (|U) + lit))/v / 2, and | J. J.), with energies E 2 = 
—iXc, E = —A, Ei = 4:Xc, E 3 = A, respectively. The energy gap \Ei — E 2 \ between the two entangled 
states is defined by &Xc ( = ^Xc)- ^ n ^ ne presence of noise, the oscillatory, approximately coherent evolution 
of the spins can be observed over several oscillation cycles provided that 2a™w c /[l + (w c |d|/c s ) 2 ] > M.. 
The energy levels will acquire effective width due to quantum noise, of order 77™ (d). 

4.3. Boson-Mediated Induced Qubit-Qubit Interaction in General Dimensions 

Let us generalize the results of the previous subsections, where wc considered the ID case with Ohmic 
dissipation within the Markovian approach. In the general case, let us consider the Markovian model and, 
again, assume that A/w c is small. We will also assume that the absolute square of the mth component of 
the spin-boson coupling, when multiplied by the density of states, can be modeled by a™ m w" m exp(— uj/w c ); 
see (14). The integration in (35) can then be carried out in closed form for any n m = 1,2,.... The induced 
interaction (43) is thus generalized to 

ti ST n ™ , ,n m 1 2 21(7i m jKe(J ±uj c \d\/c s ] 

-Hint = - 2^ a n m ^ c m <T m (J m jr- ^ 2 -|» m ■ (44) 

m— x,y ,z 

With the appropriate choice of parameters, the result for the induced interaction, but not for the noise, 
coincides with the expression for i/j nt obtained within the short-time model. The effective interaction has 
the large-distance asymptotic behavior |d| " m , for even n m , and |d| " m , for odd n m . 

In higher dimensions the structure of <?™£ in the k-space becomes important. Provided o>k,£ is nearly 

2tt tt 2 

isotropic, the integrals in (35)-(38) will include (in 3D) a factor J Q dip J Q dOsind g™ e ^ cos (k |d| cos#), 
which can be written as [/{"(w, fc|d|) - / 2 m (w, k\d\)d/d\d\] cos(fc|d|), see Eqs. (B1)-(B4) of Ref. 2 for details. 
When the dependence of /™, f™ on fe|d| is negligible, the interaction is simply H int — > H- m t\{n m }^a — 
(9/9|d|)-ffint|{n m }-»6; where a and b are sets of three integers representing the w-dependence of /"(w, A;|d|) 
and / ™( w ) /c|d|)- Otherwise, a more complicated dependence on |d| is expected. The noise superoperator 
can be treated similarly. 



2T(n 


m )Re (1 + itu c d 


/c s ) n - 




l + (uj c \d\/ Cs ) 2 


Tim 



4.4. Induced Interaction vs. Noise in the Semiconductor Impurity Qubit Model 

Let us now proceed to an example of the semiconductor qubit model in the bulk material, which gives 
usual 3D geometry for phonon propagation. We consider spins of P-impurity donor electrons as qubits. 
Dilute P-impurities are embedded in the matrix of a Si or Ge crystal kept at a sufficiently low temperature 
such that the outer electrons remain bound. The two impurities placed next to each other at distances 
of about 10-30 nm would constitute a two qubit register. The spins would then be interacting indirectly 
via the spin-orbit coupling and lattice deformation (phonons) as well as directly via the dipolc-dipolc 
electromagnetic coupling. In this section we will estimate and compare both. We will also present the 
dynamics due to the coherent vs. noise induced features and their interplay. 

Let us consider for simplicity only the LA phonon branch, £ — > k/ |k| , and assume an isotropic dispersion 
Wk,£ = c s |k|. The expression for the coupling constants is then 



One can show 2 that the cross terms, with m 7^ n, of the correlation functions Trg [A^A* (t)pB\ , depend 
on the combination g^p (g£ *)*, which is always an odd function of one of the projections of the wave 
vector. Thus, all the non-diagonal terms vanish. 

Integrating (35) and (36) with (45), one can demonstrate 2 that decoherence is dominated by the 
individual noise terms for each spin, with the typical amplitude 



27T 2 b^_ 

"TH" (1 + P) A """" 2fc B T 



V x c ' v (0) = Cb — - ; —4 coth^; , (46) 



where b = AaB/c s and Cb = B 2 /i|if 2 / ( 1 6ir 3 pha"^ c 2 ) . The interaction amplitude xJP(d) and, therefore, 
the induced spin-spin interaction, have inverse-square power-law asymptotic form for the x and y spin 
components, with a superimposed oscillation, and inverse-fifth-power-law decay for the z spin components, 



H im = J2- 2 ^(dV m a 2 m (47) 



-»i Ajin- 2b s[nbr /-i-2 , 1 „2% , 00. 2n. 2 - 1- 2 



. 2rl l + % 2 -% 4 + b 6 2 ' 

~ 4 " Cb 120(1 + vr + ° {r \ 



(T.CT, 



Here r = |d|/as and Ca = A 2 Cb/B 2 . At small distances the interaction is regular and the amplitudes 
converge to constant values, see Figure 3. The complete expressions for x™(d) and ?7™(d) are also avail- 
able. 2 

In Figure 3, we plot the amplitudes of the induced spin-spin interaction, which has the asymptotic 
behavior (47) , and noise for different values of the spin-spin separation and 6, for electron impurity spins in 
3D Si-Ge type structures. The value of b can be controlled via the applied magnetic field, b = p^H z g* a~g,/ c s . 
The temperature dependence of the noise is insignificant provided 2k^T / A -c 1. 

As mentioned earlier, the obtained interaction (47) is always accompanied by noise coming from the 
same source, as well as by the direct interactions of the spins. When the electron wave functions overlap 
is negligible, the dominant direct interaction will be the electromagnetic dipolc-dipole one, 

H EM (d) = — — . (48) 




Figure 3. The magnitude of the induced spin-spin interaction for a 3D Si-Ge type structure: the dominant 
interaction amplitude, which is the same for the x and y spin components, is shown. The arrow indicates increasing 
b values for the curves shown, with b — Aob/c s = 0.03, 0.09, 0.15, 0.21, 0.27, 0.33. The inset estimates the level 
of the noise (for 2kBTas/cs — 0.01): the bottom curve is r?™(0), with m = x,y. For b < 0.4, the amplitude ?7™(d) 
can be comparable, and its values, calculated numerically for < r < 50, are shown as long as they exceed r?™(0), 
with the top curve corresponding to the maximum value, at r = 0. 



The comparison of the two interactions and noise is shown in Figure 4. We plot the magnitude of the 
effective induced interaction (47), the electromagnetic interaction (48), measured by Hem = ho^b/^I^I 3 ^ 
and a measure of the level of noise, for P-donor electron spins in Ge. In the plot the coupling constants 
for P-donors in Ge were taken as Cb = 1.3 x 10 7 s _1 and Ca « 0, see Sec. 2. It transpires that the 
induced interaction can be considerable as compared to the electromagnetic spin-spin coupling. However 
the overall coherence-to-noise ratio is quite poor for Ge. In Si, the level of noise is lower as compared 
to the induced interaction. It is dominated primarily by the adiabatic (c 1,2 ) term. However, the overall 
amplitude of the induced terms compares less favorably with the electromagnetic coupling. 

Let us estimate the development of entanglement for the P-in-Si case, considering only the bath- 
induced effects. Taking |++) as the initial state, one expects entanglement to develop due to the WintC^c 2 
interaction term. The corresponding coupling constant is Cb = 7800 s" 1 . At distances of about 12ob this 

part of the interaction can be already well approximated by 7i\ n t ~ 3847t 2 Ca j^jp- ~ 100 s _1 . The x and y 
components of the interaction are smaller by a factor 10 -3 and will not be considered. The pure state would 

then be developing as \tp(t)) = |++) cosHint* + | )isinHi nt i. At times i# = 7r/47ii nt , Sn/AHint, ■ . ., 

the maximally entangled (Bell) states are obtained, with the first occurring at ~ 0.8 • 10~ 2 s, and then at 
intervals of ~ 1.6 ■ 10 -2 s. The actual entanglement is lowered by (i) the short-time (adiabatic) decohercncc 
term, (ii) the large-time relaxation amplitudes (rates). The latter are due to the x,y components of the 
qubit-bath coupling and are of order ~ 1 s _1 . As a result the first peak of the concurrence is lowered by a 
factor of ~ exp(— 0.8 • 10 -2 ), and the subsequent peaks decrease by ~ exp(— 1.6 • 10~ 2 ) per each cycle. 



5. DISCUSSION OF OPEN PROBLEMS 



We have studied the induced indirect exchange interaction due to a bosonic bath which also introduces 
quantum noise. At certain time scales the induced two-qubit interaction is sufficiently strong to produce 
significant coherent effects. This interaction can be a factor to be considered in designs of solid state (as well 




Figure 4. The magnitudes, measured in units of s~ , of the induced spin-spin interaction, the dipole-dipole coupling 
strength, and the level of noise for P-impurity electron spins in Ge. Here H z — 3 x 10 4 G, and low temperature, 
2fceT/A <C 1, was assumed. 



as ion-trap based) qubit registers. Even more importantly, the fact that noise-inducing environment can 
also yield coherent features in the dynamics of open quantum systems poses several interesting fundamental 
challenges. 

The usual approach to thermalization 31, 33,40-42 has been to assume that, for large enough times, the 
time evolution of the system plus bath is not just covered by the combined Hamiltonian, but is supple- 
mented by the instantaneous loss-of-memory (Markovian) approximation, which introduces irreversibility 
and imposes the bath temperature on the reduced system dynamics in the infinite-time limit, which is 
then approached as the density matrix elements assume their thermal values, according to 

PS® - e- Hs ' kT /Tt (e- H ^ kT ) , (49) 

with exponential relaxation (diagonal) and decoherence (off-diagonal) rates defining the time scales Ti,2- 
However, it turns out that the traditional phenomenological no-memory approximations, yielding ther- 
malization, the Fermi golden rule for the transition rates, etc., assume in a way too strong a memory 
loss: they erase the possible bath contributions to the coherent part of the dynamics at shorter times, 
such as the Lamb shift for a single system as well as the induced RKKY interactions for a bi-partitc 
system. Indeed, while relaxation leading to (49) is driven by the "on-shell" exchanges with the bath, it 
is the memory of (correlation, entanglement with) the bath modes that drives, via virtual exchanges, the 
induced interaction. Actually, the "on-shell" condition, imposed by the so-called secular approximation, 7 
also underestimates additional decoherence at short times — the "pure" or "adiabatic" contribution to the 
off-diagonal dephasing — that has thus been estimated by using other approaches. 28, 29,43-45 

Perhaps the simplest way to recognize the ambiguity is to ask if the Hamiltonian in (49) should 
have included the bath- induced interaction terms (not shown)? Should the final Boltzmann distribution 
correspond to the energy levels/basis states of the original "bare" system or the one with the RKKY- 
intcraction/Lamb shifts, and more generally, a bath-renormalized, "dressed" system? 

There is presently no consistent treatment that will address in a satisfactory way all the expected 
physics of the bath-mode effects on the dynamics. The issue is partly technical, because we are after a 
tractable, rather than just a "foundational" answer. It is well accepted that the emergence of irreversibility 



cannot be fully treated within tractable and calculationally convenient approaches derived directly from the 
microscopic dynamics: phenomenology has to be appealed to. However, even allowing for phcnomenological 
solutions, most of the known tractable Markovian-approximation-involving schemes that allow for the 
emergence of the indirect exchange interaction, treat thermalization in a cavalier way, yielding typically 
the noise term corresponding to T-dependent relaxation, but in the strict t — * oo limit resulting in the 
completely random (T = oo) density matrix (proportional to the unit operator). This then avoids the 
issue of which H$ should enter in (49). And, as mentioned, the established schemes that yield a more 
realistic, thcrmalizcd density matrix at t = oo, lose some intermediate- and short-time dynamical effects. 

Thus, we have discussed the challenges in formulating unified treatments that will cover all the (or just 
most of the interesting) dynamical effects, over several time scales, from short to intermediate times (for 
induced interaction effects and pure decoherence) to large times (for the onset of thermalization), while 
providing a tractable calculational (usually perturbativc, many-body) scheme. This discussion also alludes 
to several other interesting conceptual challenges in the theory of open quantum systems. 

Let us presently comment on the issue of the bath-mode interactions with each other, as well as with 
impurities, the latter particularly important and experimentally relevant 46 ' 47 for conduction electrons as 
carriers of the indirect-exchange interaction. Indeed, the traditional treatment of open quantum systems 
has assumed noninteracting bath modes. When the bath mode interactions had to be accounted for, the 
added effects were treated cither pcrturbativcly, 48 or, for strong interaction, such as Luttinger-liquid elec- 
trons in a ID channel, the collective excitations were taken 49 as the new "bath modes." Generally, however, 
especially in Markovian schemes, one has to seek approaches that do not involve certain double-counting. 
Indeed, the assumption that the bath modes are at a fixed temperature, could be possibly considered 
as partially accounting for the effects of the mode- mode and mode- impurity interactions, because these 
interactions can contribute to thermalization of the bath, on par with other influences external to the bath. 
This may be particularly relevant for phonons that always have strong anharmonicity for any real material. 
In a way this problem fits with the previous one: we are dealing with effects that can be, on one hand, 
modeled by added terms in the total Hamiltonian but on the other hand, may be also mixed in the process 
of thermalization that is modeled by actually departing from the Hamiltonian description and replacing 
it with Liouville equations that include noise effects. While all this sounds somewhat "foundational," 
recent experimental advances, interestingly, bring these challenges to the level of application that requires 
tractable, albeit perhaps phcnomenological solutions that can be directly confronted with experimental 
data. 

There are other interesting topics to be considered, for instance, the question of whether additional 
sources of quantum noise are possible? For instance, it has been recently established 50 that potential 
difference between two leads (reservoirs, or baths, of electrons) can be a source of quantum noise with the 
potential difference playing the role of the temperature parameter. 

As an example of a more practical issue to be investigated, let us mention the possible effect of the sam- 
ple geometry on phonon and conduction electron induced relaxation and interactions. The one-dimensional 
aspect of the electron gas in a channel has already been explored. 49 Indeed, electrons are easy to confine 
by gate potentials. The situation for phonons, however, is less clear: can geometrical effects modify, and 
particularly reduce, their quantum-noise generation capacity, or change the induced interactions? Our 
preliminary studies reviewed in this article, seem to indicate strong overall geometry dependence of the 
exchange interactions. However, the situation is not entirely clear, especially for the strength of the noise 
effects, and requires a full scale exploration because recent experiments with double dot nanostructures 
in Si membranes 51 suggest that true nanosize confinement (due to the sample dimensions) of otherwise 
long-wavelength modes (in the transverse sample dimensions) is now possible and will have dramatic ef- 
fect on the phonon spectrum and, as a result, on those physical phenomena that depend on the phonon 
interactions with electron spins. 
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